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ABSTRACT 


The  combination  of  cartographic  and  hyposographic  data  with  radar 
data  can  generate  the  location  of  terrain  features  in  the  ground 
reference  system.  The  significance  of  selected  data  collection  parameters 
used  in  this  task  is  investigated  in  three  phases: 

1)  The  development  of  deterministic  radar  models  from  given 
digital  terrain  models  and  varying  radar  parameters. 

2)  The  perturbation  of  the  terrain  feature  coordinates  which  result 
from  errors  in  radar  position,  radar  orientation,  range,  and 
resolution. 

3)  A  factorial  analysis  of  the  terrain  feature  errors  to  establish 
the  significance  of  the  main  fixed  factors  and  their  interactions. 
The  main  factors  are  azimuth  from  radar  to  terrain  feature,  range, 
resolution,  terrain  height  and  coordinate. 

The  results  indicate  that  the  order  of  greatest  to  least  significance 
for  the  factors  studied  is  coordinate,  range,  resolution,  terrain  height 
and  azimuth.  Additional  results  are  given  in  fourteen  figures  and  eleven 
tables.  An  appendix  is  included  on  the  projection  equations  for  real- 
aperture  radar. 


i  i  i 


PriECEOlNG  PACK  bunk -not  Hi>£D 


THE  SIGNIFICANCE  OF  DATA  COLLECTION  PARAMETERS  WHEN 
UTILIZING  RADAR  AND  DIGITAL  TERRAIN  MODELS 
TO  LOCATE  TERRAIN  FEATURES 


1.0  INTRODUCTION 
1.1  BACKGROUND 

The  continuing  improvement  in  radar  systems  coupled  with  their 
unique  capabilities  of  information  collection  indicate  that  they  will 
play  an  increasing  role  in  future  mapping  systems.  The  ultimate  role 
or  path  to  that  role  has  yet  to  be  determined. 

Considerable  experimentation  has  taken  place  during  the  past  decade 
with  the  goal  of  determining  the  value  of  radar  as  a  mapping  or  terrain 
feature  positioning  tool.  The  geometry  applicable  to  radar  mapping 
has  been  well  documented,  and  the  mathematical  models  have  been  developed 
with  the  thought  of  combining  the  radar  sensor  with  other  types  of  remote 
sensors.  The  most  significant  combinations  of  sensors  have  yet  to  be 
determined. 

The  digital  terrain  model  (DTM)  is  an  emerging  information  source 
which  will  have  an  impact  on  the  use  of  radar  as  a  mapping  tool.  Ground 
truth  is  an  important  part  of  any  mapping  task  and  the  availability  of 
DTM  data  can  be  used  in  many  ways  to  provide  a  geometric  framework  for 
the  radar  data.  The  optimum  combinations  of  these  two  data  sources  will 
not  be  easy  to  ascertain  and  this  topic  will  need  considerable  study. 
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1.2  CONCEPT  AND  PURPOSE  OF  THE  INVESTIGATION 


The  overall  concept  of  the  mapping  task  considered  in  this  study 
is  to  combine  the  capability  of  the  radar  with  the  DTM  to  produce  the 
location  of  terrain  features  in  the  X,  Y,  and  Z  DTM  coordinate  system. 

The  assumption  is  that  a  DTM  is  available  for  the  area  of  Interest.  A 
radar  sensor  is  to  be  flown  in  the  vicinity  of  the  area  of  concern  in 
order  to  collect  imagery. 

Assuming  a  given  sensor,  whose  general  capability  is  known,  the 
object  is  to  determine  the  most  efficient  operational  use  of  the  system. 
The  determination  of  target  positions  is  a  function  of  several  para¬ 
meters  whose  values  are  known  within  certain  limits.  If  the  values  of 
the  data  collection  parameters  are  changed,  will  the  results  also  be 
changed?  Some  of  the  parameters  can  be  varied  mathematical ly  for  strong¬ 
er  functional  relationships,  however,  some  depend  solely  on  the  original 
data  collection  procedures.  Accordingly,  the  purpose  of  this  study  was 
to  select  several  important  data  collection  parameters  and  to  study 
their  significance  in  the  location  of  terrain  features  when  radar  and 
DTM  data  are  combined.  The  determination  of  the  relative  importance  of 
data  collection  parameters  should  aid  in  the  design  of  successful  data 
collection  efforts.  For  example,  Figure  1  schematically  demonstrates 
the  effect  one  can  expect  from  the  variation  of  range,  azimuth,  and 
terrain  height  during  the  data  collection  process.  The  standard  devia¬ 
tions  are  for  the  conditions  imposed  in  this  study. 


2 


The  general  procedure  employed  was  to  develop  deterministic  models, 
to  perturb  the  selected  parameters  in  a  random  fashion,  to  establish 
errors  for  specific  parameter  combinations  and  to  evaluate  the  results 
with  a  factorial  analysis.  Section  2  elaborates  on  the  details  of  the 
procedure  and  section  3  states  and  discusses  the  results.  Section  4 
gives  conclusions  and  recommendations.  Appendices  A  and  B  provide 
details  on  the  computational  and  statistical  procedures. 


S 
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R-Radar 


p -Terrain  Feature 


2.0  EXPERIMENTAL  PROCEDURE 
2.1  GENERAL 

This  section  outlines  the  procedure  employed  during  the  study.  It 
consists  of  several  sub-sections,  namely;  the  definition  of  the  deter¬ 
ministic  model  which  served  as  the  basis  of  the  study,  the  error  pro¬ 
pagation  phase  of  the  effort  which  perturbed  the  deterministic  model 
in  a  random  fashion,  and  the  factorial  analysis  of  the  errors  introduced 
by  perturbing  the  selected  parameters. 


2.2  DETERMINISTIC  MODEL 


Two  coordinate  systems  are  utilized  in  setting  up  the  model.  The 
first  is  the  ground  coordinate  reference  system  which  is  the  basis  for 
the  DTM.  In  this  instance  the  RADC  DME  Navigation  Range  in  the  vicinity 
of  Rome,  N.Y.  was  selected  for  the  study.  Cartographic  and  hypsographic 
data  have  been  collected  for  this  area  and  a  DTM  exists.  A  local  radar 
coordinate  system  also  has  to  be  established  and  an  expression  adopted 
for  the  transformation  of  one  system  into  the  second  and  vice  versa. 
Figure  2  shows  an  example  of  the  relative  position  of  the  two  reference 
systems.  Equation  1  adopted  from  Franz  Leberl  states  the  mathematical 
model  which  is  the  basis  for  the  study.  (Leberl,  1976). 
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A  =  Orientation  matrix  in  terms  of  k,$,w 
as  normally  adopted  in  photogrammetry. 

r  =  Range  from  radar  to  DTM  point  P. 

t  =  Squint  angle 

fi  =  Elevation  angle 

A  total  of  27  unique  radar  scenes  were  generated  with  equations  (1) 
and  (2).  X,  V,  and  Z  coordinates  for  twenty-one  points  were  generated  in 
each  scene  which  nominally  covered  an  area  in  the  DTM  of  5  km,  by  5  km. 
This  experiment  was  repeated  three  times  for  each  of  the  27  scenes.  Fig¬ 
ure  3  shows  the  positions  of  the  selected  21  control  points  in  the  DTM. 

Recalling  that  the  purpose  of  the  study  is  to  determine  the  relative 
importance  of  radar  data  collection  parameters,  the  next  task  is  to 
select  the  main  parameters  or  factors  for  study.  The  parameters  studied 
in  this  effort  were  range,  azimuth  cr  radar  direction  which  was  angle 
kappa  in  this  application,  the  resolution  of  the  radar,  the  height  of  the 
DTM  point,  and  the  particular  coordinate  for  the  point  of  interest. 

When  locating  random  points  in  a  racar  scene  and  transforming  the  ,'adar 
location  into  a  meaningful  position  in  the  DTM,  any  or  all  of  these  para¬ 
meters  may  be  important.  The  task  ii  to  determine  how  the  relative  val¬ 
ues  of  these  parameters  affect  the  accuracy  of  the  final  point  position 
in  the  DTM.  Are  their  particular  combinations  to  either  avoid  or  encour¬ 
age  for  the  sake  of  accurate  mapping  or  position  determination? 

The  total  experiment  is  summarized  in  Figure  4.  For  each  of  the 
computations  noted  in  Figure  4,  the  values  of  range  (r),  squint  angle 
(t)  and  elevation  angle  (n)  were  computed  when  the  values  Xp,  Y  ,  Z  , 

Xp,  Yr,  ZR,  k,  4> ,  and  w  were  provided. 


7 


8 


Formulas  used  for  the  computation  follow: 


RANGE  (r)  =  [(Xp  -  XR)2  +  (Yp  -  YR)2  +  (Zp  -  ZR)2  ** 

SQUINT  ANGLE  (i)  =  SIN-1  \*U  (X  -XR)  +  A12  (Yp-YR)  +  A31  (Zp-Z  ) 

I  - ' - I 

1—  r  -i 

ELEVATION  ANGLE  (ft)  =C0S_1|  ^  +  (Yp'V  +  ^  ^W"] 

L  ^  j 


2.3  ERROR  MODELS 

The  computed  DTM  coordinates  (X„  1)  are  a  function  of  XD,VD,ZD 

P  >  P  5  P  KKK, 

k,<t>,w,r,T,  and  ft.  Under  operational  conditions  these  values  must  be 
provided  either  directly  or  indirectly.  Each  value  will  have,  at  least, 
a  random  error  associated  with  it.  Previous  studies  and  experimental 


measurement  in  this  area  by  the  authors  and  others  provide  us  with  rea¬ 
sonable  estimates  for  the  appropriate  random  error  to  be  applied  for 
each  parameter.  (Brock,  Zulch  1979). 

The  position  of  the  radar  is  provided  by  DME  navigation.  It  s 


estimated  that  XR=±9  m.,  YR=±3  m.,  anu  ZR-±20  m. 


Stabalizeo  platfvrm., 


are  estimated  to  provide  k,<j>,  ana  w  to  -0.1°.  The  error  in  range  . s  set 


at  ±5  m.  One  of  the  additional  factors  was  radar  resolution.  Two  cases 


were  assumed.  The  first  assumed  that  the  minimum  element  of  information 


had  a  diameter  of  5  m.  and  the  second  case  assumed  that  the  minimum 


element  size  was  10  m.  For  a  range  of  200  km,  for  instance,  this  would 
then  place  the  random  error  for  x  and  u  at  i0.0015°  for  the  5  m.  case 
and  ±0.003°  for  the  10  m.  case. 
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Utilizing  the  above  standard  deviations  and  a  random  number  gener¬ 
ator,  random  errors  were  applied  to  each  of  the  9  parameters  utilized  in 
the  computation  of  the  Xp,Yp,  and  zp  coordinates.  This  provided  per¬ 
turbed  values  for  X  ,Y  ,  and  Z  .  The  perturbed  values  were  subtracted 

r  r  r 

from  the  true  values  previously  assigned.  The  resulting  errors  provided 
the  basis  for  the  factorial  analysis  portion  of  the  study.  See  Table  1 
for  an  example  of  true  and  perturbed  values. 

2.4  FACTORIAL  ANALYSIS 

The  main  objective  of  the  effort  was  to  demonstrate  the  importance 
of  azimuth,  range,  resolution,  terrain  height  and  coordinates  in  the 
utilization  of  radar  and  DTM  data  to  provide  DTM  positions  of  random 
points  of  interest. 

A  factorial  analysis  was  selected  since  it  not  only  allows  one  to 
examine  the  effects  of  each  parameter  but  the  effects  of  the  interac¬ 
tions  of  the  factors  as  well.  In  addition  all  values  may  be  used  to 
investigate  each  of  the  factors.  The  possibility  of  ascribing  variation 
to  the  wrong  factor  is  reduced  since  all  factors  are  contained  in  the 
model.  A  summary  of  the  main  experimental  factors  and  the  levels  of  each 
is  given  in  Figure  4. 
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TABLE  1  -  SAMPLE  OF  TRUE  AND  PERTURBED  VALUES  FOR 

PARAMETERS  AND  A  TERRAIN  FEATURE  OF  INTEREST 

PARAMETER  TRUE  VALUE  PERTURBED  VALUE 


XR 

302,000 

M. 

302,010 

M. 

yr 

523,000 

M. 

522,997 

M. 

ZR 

12,192 

M. 

12,187 

M. 

k 

1.570796 

Rad. 

1.569871 

Rad 

<f> 

0.0  Rad. 

0.002028 

Rad 

w 

0.0  Rad 

0.001062 

Rad 

r 

198,944 

M. 

198,936 

M. 

T 

0.000000 

Rad. 

-0.000024 

Rad 

a 

1.512017 

Rad. 

1.512037 

Rad 

TERRAIN 

FEATURE 

V 

500,600 

M. 

500,590 

M. 

¥p 

523,000 

M. 

523,162 

M. 

z„ 

505 

M. 

2':)  4 

M. 
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3.0  RESULTS 


3.1  GENERAL 

The  results  will  be  presented  in  three  parts.  First,  sample  results 
will  be  given  for  the  generation  of  deterministic  models.  Second,  the 
results  of  several  error  models  will  be  listed  and  last  the  factorial 
analysis  results  will  be  stated.  Where  appropriate,  graphical  results  will 
also  be  given.  The  last  portion  of  this  section  will  deal  with  a  discus¬ 
sion  of  the  results  as  found  for  the  conditions  of  this  study. 

3.2  DETERMINISTIC  MODEL  RESULTS 

As  previously  stated  a  total  of  27  unique  radar  scenes  were  mathemat¬ 
ically  created.  First,  the  radar  was  directed  along  three  different  azi¬ 
muths  with  each  azimuth  being  used  for  three  distinct  ranges.  Figure  5 
shows  the  nine  radar  positions  created  by  these  two  variables.  For  each 
of  the  nine  positions,  scenes  having  three  different  average  heights  were 
created  thus  making  27  scenes.  A  set  of  21  ficticious  rays  were  created 
for  each  of  the  27  scenes.  The  range,  squint  angle  and  elevation  angle 
were  calculated  for  each  ray.  Figure  6  shows  partial  results  for  one  of 
the  27  scenes. 

3.  3  RESULTS  OF  ERROR  MODELS 

As  discussed  in  Section  2.3  the  21  ficticious  rays  created  for  each 
scene  were  randomly  perturbed  by  introducing  random  errors  into  the  radar 
position  and  orientation  and  the  ray  distance  and  direction.  For  a  given 
set  of  parameters  a  set  of  X,  Y,  and  Z  error  components  were  calculated. 
Tables  2  through  6  are  five  sets  of  errors  selected  from  the  162  sets 
which  were  computed.  These  error  sets  provide  examples  of  noteworthy 
results  which  are  discussed  in  Section  3.5. 
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RADAR 


FIGURL  5.  RADAR  POSITIONS  USED  IN  THE  STUDY 


DTM 

POINT 

SQUINT 

ANGLE  (t)° 

ELEVATION 
ANGLE  (0)0 

5 

102,138 

-0.79 

83.15 

6 

99,332 

-0.82 

82.95 

8 

102,138 

0.79 

83.15 

9 

99,331 

0.82 

82.95 

FIGURE  6.  PARTIAL  RESULTS  FOR  RADAR  POSITION  8  AS 


SHOWN  IN  FIGURE  5. 


Table  2  -  PROPAGATED  ERRORS  FOR  CONTROL  POINT 

COORDINATES-M. 

Azimuth  - 

45°,  Range  200  km., 

Resolution 

-5  m. 

Terrain  Height  - 

0  m. 

Control 

Point 

Coordinates 

Resultant 

No, 

X 

Y 

Z 

(x2  +  y2  f  z2) 

1 

251.86 

-361.70 

-963.87 

1059.89 

2 

-399.83 

324.90 

-480.73 

704.64 

3 

-611.52 

658.00 

782.70 

1191.44 

4 

-458.83 

530.01 

391.82 

803.08 

5 

468.66 

-501.44 

-772.95 

1033.70 

6 

25.09 

-  12.61 

534.84 

535.58 

7 

-238.15 

276.26 

-200.15 

416.05 

8 

563.18 

-509.07 

828.64 

1123.82 

9 

-473.49 

463.31 

-417.57 

783.08 

10 

488.53 

-356.36 

974.50 

1146.87 

11 

-500.96 

477.97 

-421.11 

810.40 

12 

-568.16 

625.99 

746.33 

1127.69 

13 

-451.65 

382.80 

-609.50 

849.71 

14 

-558.00 

583.47 

-  23.99 

807.70 

15 

707.39 

-684.52 

135.20 

993.60 

16 

2.05 

-  37.60 

-459.08 

460.62 

17 

193.21 

-172.96 

421.29 

494.70 

18 

-148.85 

205.69 

696.21 

741.06 

19 

-164.81 

101.04 

-796.52 

819.64 

20 

695.21 

-753.83 

-354.92 

1085.15 

21 

689.66 

-587.81 

896.70 

1274.84 

Average 

-23.30 

31.03 

43.23 

58.09 

Std.  Deviation 

462.32 

459.41 

o24.25 

902.49 
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Table  3  -  PROPAGATED  ERRORS  FOR  CONTROL  POINT  COORDINATES-M. 
Azimuth  -  45°,  Range  -  50  km.,  Resolution  -5m. 

Terrain  Height  -  0  m. 


Control  Point 

Coordinates 

Resultant 

No. 

X 

Y 

Z 

(X2  +  Y^  +  Z2) 

I 

61.68 

-  2.74 

227.35 

235.58 

2 

41.57 

-  83.88 

-104.46 

140.27 

3 

78.54 

-147.57 

-181.23 

246.56 

4 

-95.74 

34.73 

-246.95 

249.38 

5 

-42.64 

36.83 

114.57 

127.67 

6 

-27.56 

-  45.13 

-173.29 

181.18 

7 

107.27 

-220.22 

-138.95 

281.62 

8 

44.21 

-  39.50 

81.89 

101.10 

9 

-73.99 

86.87 

5.29 

114.23 

10 

77.91 

10.71 

184.49 

200.55 

11 

-110.47 

42.34 

-105.30 

158.38 

12 

128.95 

-  92.98 

197.62 

253.63 

13 

-  40.40 

129.48 

192.84 

235.76 

14 

-  14.27 

-  67.76 

-186.96 

199.37 

15 

-109.22 

42.54 

-260.00 

285.20 

16 

-  37.55 

54.78 

-  43.93 

79.63 

17 

19.78 

-  9.09 

70.77 

74.04 

18 

-  3.76 

15.15 

125.09 

126.06 

19 

134.49 

-  70.51 

205.62 

255.61 

20 

-  36.62 

21.63 

-133.00 

139.63 

21 

-  65.05 

87.16 

64.40 

126.40 

Average 

1.77 

-  10.34 

-  4.91 

11.58 

Std.  Deviation 

74.47 

81.31 

159.76 

194.11 
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Table  4  -  PROPAGATED  ERRORS  FOR  CONTROL  POINT  COORDINATES-M. 
Azimuth  -  90°,  Range  -  50  km.,  Resolution  -5m. 


Terrain  Height  -  0  m. 

Control  Point 

Coordinates 

..Resultant 

No. 

X 

Y  Z 

(X2  +  Y*  +  Z2) 

1 

75.67 

-  19.71 

272.29 

283.30 

2 

-  26.76 

-161.52 

-  79.13 

181.84 

3 

53.82 

186.73 

35.14 

197.48 

4 

-  39.09 

108.55 

-190.51 

222.72 

5 

-  11.04 

-251.48 

29.42 

253.44 

6 

-  10.97 

96.68 

-110.03 

146.88 

7 

16.10 

86.07 

100.21 

133.08 

8 

82.05 

-  75.98 

241.17 

265.83 

9 

-  51.22 

24.95 

-147.83 

158.43 

10 

-  14.10 

98.65 

70.22 

121.91 

11 

-  1.39 

-  61.22 

-  40.72 

73.54 

12 

23.88 

84.53 

143.57 

168.31 

13 

31.51 

36.44 

229.79 

234.79 

14 

-  59.16 

30.44 

-178.75 

190.73 

15 

31.29 

170.94 

117.80 

209.94 

16 

-  55.89 

181.36 

-190.76 

269. OS 

17 

41.50 

159.50 

273.09 

318.97 

18 

24.49 

-  29.48 

208.80 

212.29 

19 

29.79 

-120.26 

146. S4 

192.12 

20 

6.75 

-195.79 

.  32.87 

198.64 

21 

-  7.36 

115.39 

-  41.03 

122.69 

Average 

6.66 

22.13 

43.93 

49.64 

Sid.  Deviation 

39.32 

124.72 

152.00 

200.51 
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Table  5  -  PROPAGATED  ERRORS  FOR  CONTROL  POINT  COORDINATES-M. 
Azimuth  -  45°,  Range  -  50  km..  Resolution  -5m. 

Terrain  Height  -  2000  m. 


Control  Point 
No. 

X 

Coordinates 

Y 

Z 

1 

175.31 

-117.24 

184.55 

280.25 

2 

-119.38 

158.96 

130.09 

237.58 

3 

9.93 

53.22 

258.75 

264.35 

4 

219.38 

-196.15 

140.38 

326.05 

5 

-  27.31 

77.12 

236.37 

250.13 

6 

-160.57 

57.03 

-274.87 

323.40 

7 

-  73.52 

47.75 

-224.67 

241.17 

8 

-  28.81 

-  19.45 

-109.71 

115.08 

9 

84.48 

-204.11 

-178.65 

284.10 

10 

-123.69 

103.50 

-  86.44 

182.98 

11 

-  15.55 

78.26 

276.94 

288.21 

12 

26.16 

18.63 

25.69 

41.13 

13 

72.70 

-102.01 

-  76.65 

146.86 

14 

58.26 

-  75.34 

-  59.56 

112.33 

15 

156.55 

-134.50 

50.44 

212.47 

16 

192.61 

-145.81 

60.39 

249.01 

17 

10.72 

-  58.44 

-248.97 

255.96 

18 

-146.65 

167.19 

75.49 

234.86 

19 

77.90 

-  46.74 

137.00 

164.38 

20 

79.29 

-155.06 

-198.27 

263.90 

21 

100.16 

-  82.92 

80.32 

152.84 

Average 

27.05 

-  27.41 

9.46 

39.65 

Std.  Deviation 

108.73 

110.25 

168.22 

228.64 
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Table  6  -  PROPAGATED  ERRORS  FOR  CONTROL  POINT  COORDINATES-M. 
Azimuth  -  45°,  Range  -  50  km..  Resolution  -  10  m. 
Terrain  Height  -  0  m. 


Control  Point 

No.  X 


Coordinates 

Y 


Resultant 

Z  (X2  +  Y2  +  Z2)’* 


1 

-132.38 

2 

-170.43 

3 

162.97 

4 

-112.00 

5 

-  28.60 

6 

6.11 

7 

29.34 

8 

-  22.00 

9 

169.34 

10 

22.43 

11 

62.65 

12 

145.87 

13 

38.87 

14 

35.68 

15 

227.89 

16 

135.80 

17 

-  6.90 

18 

67.37 

19 

-198.92 

20 

84.06 

21 

-  32.08 

Average 

23. 

,10 

Std.  Deviation 

110, 

.51 

190.30 

154.91 

-131.66 

115.17 

-  7.06 

-  7.05 

-  42.27 
14.56 

-151.47 

-103.11 

-  96.92 
-115.61 

41.58 

11.53 

-178.53 

-  64.43 

-  58.08 

-  96.69 
151.72 

-145.40 

34.52 


-  23.04 
1U5.57 


256. 

18 

345. 

49 

-115. 

42 

257. 

61 

120. 

.20 

241. 

54 

-  55. 

15 

169. 

85 

-161. 

,37 

164. 

04 

20. 

,43 

22. 

46 

-  82. 

,77 

97. 

,46 

81. 

.51 

85. 

,67 

48. 

,74 

232. 

,37 

-225. 

.66 

249. 

,11 

-271 

.02 

294. 

.57 

140. 

.22 

233. 

.04 

224, 

.91 

232. 

.00 

137, 

.85 

142 

.86 

52, 

.40 

294. 

.20 

186, 

.40 

239, 

.45 

-123 

.94 

137 

.05 

-  32, 

.70 

122, 

.22 

-222 

.39 

334 

.73 

-166 

.32 

237 

.78 

-122 

.15 

130 

.93 

-  14. 

153.43 

216.56 
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The  computations  in  Table  3,  for  example,  were  repeated  three  times, 
each  with  a  different  set  of  random  errors.  These  three  computations 
were  referred  to  as  experiments.  Table  7  shows  the  average  error  results 
of  these  three  computations  for  the  parameter  values  given  in  Table  3. 


TABLE  7.  AVERAGE  ERROR  RESULTS  FOR  THREE 

EXPERIMENTS  (Table  3  parameters  used) 


Experiment 

X 

Coordinate  -  m. 

Y 

Z 

1 

1.76 

-10.34 

-  4.92 

2 

44.87 

-43.37 

23.36 

3 

-46.88 

43.40 

-  7.33 

Average 

-  0.08 

-  3.44 

3.70 

Std.  Deviation 

45.90 

43.79 

17.07 

This  is  an 

important  point  of 

the  study  since  the 

test  for  signifi- 

cance  is  against 

the  experimental 

error.  If  one  cannot 

repeat  a  coordi- 

nate  determination  to  better  than  ±45  m. ,  then  a  parameter  cannot  be 
judged  to  be  significant  unless  the  error  produced  by  the  variation  of 
the  parameter  is  statistically  larger  than  ±45  m. 

3.4  RESULTS  OF  ANALYSIS  OF  VARIANCE 

The  basic  null  hypothesis  (H0)  for  the  analysis  of  variance  (ANOVA) 
model  in  this  study  is  one  of  no  significance  between  the  factors  or 
interaction  of  the  factors  at  the  10%  level.  Thus,  if  the  analysis  shows 
significance,  it  is  indicating  that  one  could  expect  a  Type  I  error  10 
times  out  of  100  in  the  long  run. 
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The  results  of  the  Factorial  Analysis  are  given  in  Table  8,  where  the 

fixed  factors  of  concern  are  as  follows: 

A  -  Azimuth  D  -  Terrain  Height 

B  -  Range  E  -  DTM  Coordinate 

C  -  Resolution 

The  F  test  for  significance  depends  upon  the  variable  introducing  a 
variance  which  is  significantly  larger  than  that  generated  by  completing 
the  experiment  independently  three  times.  It  is  seen  that  only  the  DTM 
coordinate  is  significant.  The  order  of  importance  for  the  main  factors 
from  most  to  least  is  DTM  coordinate,  range,  resolution,  terrain  height 
and  azimuth. 

The  ANOVA  computation  also  provides  the  average  errors  for  all  lev¬ 
els  of  the  five  main  factors.  See  Table  9  for  these  values. 

Referring  to  the  coordinate  results  it  is  seen  that  the  significance 
indication  can  be  attributed  to  the  large  variations  in  the  V  coordinate. 

Considering  the  interaction  of  any  two  of  the  factors  it  is  seen 
that  the  combination  of  range  and  terrain  height  provide  significance. 

Most  of  the  variance  can  be  traced  to  the  case  of  200  KM  range  and  20(0  v. 
terrain  height.  (Table  10).  A  look  at  throe  factor  interactions  unoc 
significance  for  range,  terrain  height  and  coordinate  combinations,  'ne 
greatest  variance  is  provided  at  the  200  KM  range  and  in  particular  the 
200  KM  range,  2000  M.  terrain  height  and  Z  coordinate.  (Table  11). 

Four  factor  interactions  have  one  significant  combination,  namely;  azimuth, 
resolution,  terrain  height  and  coordinate.  The  intermediate  level  of  10 
M.  resolution,  500  M.  terrain  height,  Z  coordinate  and  90°  azimuth  gives 
the  greatest  variance. 
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TABLE  8  -  ANALYSIS  OF  VARIANCE  RESULTS 


Source 

Degrees 

Of  Freedom 

Sum  of  Squares 

(m.)2 

Mean  Square 
(m.)2 

F.  Test 

A  -  Azimuth 

2 

6702 

3551 

0.59 

B  -  Range 

2 

23477 

11738 

2.06 

C  -  Resolution 

1 

6342 

6343 

1.12 

D  -  Terrain  Ht. 

2 

11168 

5584 

0.98 

E  -  Coordinate 

2 

29162 

14581 

2.56* 

AB 

4 

5618 

1405 

0.25 

AC 

2 

7361 

3680 

0.65 

AD 

4 

23147 

5787 

1.02 

BC 

2 

2741 

1371 

0.24 

BD 

4 

48529 

12132 

2.13* 

CD 

2 

1452 

726 

0.13 

AE 

4 

23263 

5816 

1.02 

DE 

4 

33176 

8294 

1.46 

BE 

4 

34947 

8737 

1.54 

CE 

2 

524 

262 

0.05 

ABC 

4 

3556 

889 

0.16 

ABD 

8 

22820 

2853 

0.50 

BCD 

4 

6423 

1606 

0.28 

ACD 

4 

23088 

5772 

1.02 

CDE 

4 

21207 

5302 

0.93 

DEA 

8 

26420 

3303 

0.58 

ABE 

8 

30044 

3755 

0.66 

BDE 

8 

86500 

10813 

1.90* 

ACE 

4 

18628 

4657 

0.82 

BCE 

4 

13116 

3279 

0.58 

ABCD 

8 

10640 

1330 

0.23 

BCDf 

8 

33907 

4238 

0.75 

COEA 

8 

85927 

10741 

1.89* 

ABDt 

16 

40329 

2521 

0.44 

ABCE 

8 

27373 

3422 

0.60 

ABCDE 

16 

92878 

5805 

1.02 

Experimental 

324 

1842264 

5686 

Error 

Total 

485 

2642728 

*  Significant  at  10%  level. 
All  others  not  significant. 
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TABLE  9  -  AVERAGE  ERRORS  IN  METERS  FOR  MAIN  FACTORS 
(MAIN  FACTORS  ARE  LISTED  IN  ORDER  OF  SIGNIFICANCE) 


Order  Factor  ( EX . ) ^*  Ave.  X- 

of  results  Values  /  ?\  . 1 

(ro*  )  (I  m. ) 


Coordinate  Results 


1 

Z 

63,360 

1.6 

2 

X 

191,784 

2.7 

3 

Y 

6,255,196 

15.4 

Range  Results 

1 

100 

km 

27,626 

1.1 

2 

50 

km 

40,453 

1.2 

3 

200 

km 

5,521,278 

14.5 

Resolution  Results 

1 

5  m.  78,144 

1.2 

2 

10  m.  4,142,361 

8.4 

Terrain 

Height  Results 

1 

0  m. 

21,734 

0.9 

2 

500  m. 

504,628 

4.4 

3 

2000  m. 

3,069,059 

10.8 

Azimuth  Results 

1 

45° 

20,021 

0.9 

2 

60° 

349,957 

3.6 

3 

90° 

2,501,950 

9.8 

*  X.  =  Average  error  in  meters  for  21  point  scene. 
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TABLE  10  -  AVERAGE  ERRORS  FOR  INTERACTIONS  BETWEEN 


Order  of 
Result 


TWO  OF  THE  MAIN  FACTORS. 

(SELECTED  COMBINATIONS  ARE  GIVEN  IN  ORDER  OF  MOST 
SIGNIFICANT  TO  LEAST  SIGNIFICANT). 


Combination 

Values 


(EX^2  Ave.  X. 

(m*2)  (t  m.) 


Range  and  Terrain  Height  Results 


Range 

Terrain  Ht. 

50  km 

500  m. 

3 

0.1 

100  km 

2000  m. 

11,137 

2.0 

50  km 

0  m. 

11,437 

2.0 

100  km 

0  m. 

14,755 

2.2 

200  km 

0  m. 

17,662 

2.5 

50  km 

2000  m. 

93,817 

5.7 

200  km 

500  m. 

101,721 

5.9 

100  km 

500  m. 

154,618 

7.3 

200  km 

2000  m.  4 

,681,606 

40.0 

Range  and  Coordinate  Results 


Range 

50  km 

Coordinate 

X 

378 

0.4 

50  km 

Y 

598 

0.5 

100  km 

X 

4,603 

1.3 

100  km 

Z 

6,850 

1.5 

50  km 

Z 

60,035 

4.5 

100  km 

Y 

100,374 

5.9 

200  km 

X 

151,727 

7.2 

200  km 

Z 

335,815 

10.7 

200  km 

Y 

4,664,588 

40.0 
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TABLE  10  -  (CONT'D) 

Terrain  Height  and  Coordinate  Results 


Terrain  Ht. 

500  m. 

Coordinate 

Z 

3,348 

1.1 

2000  m. 

X 

36,290 

3.5 

500  m. 

Y 

81,175 

5.3 

500  m. 

X 

135,128 

6.8 

0  m. 

X 

378,261 

11.4 

0  m. 

Z 

561,545 

13.9 

2000  m. 

Z 

889,652 

17.5 

2000  m. 

Y 

998,319 

18.5 

0  m. 

Y 

1,481,013 

22.5 

Azimuth  and  Terrain  Height 

Azimuth 

Terrain  Height 

45° 

0  m. 

5080 

1.3 

45° 

500  m. 

48,293 

4.1 

60° 

500  m. 

63,594 

4.7 

60° 

0  m. 

83,260 

5.3 

90° 

0  m. 

133,005 

6.8 

45° 

2000  m. 

187,077 

8.0 

60° 

2000  m. 

308,249 

10.3 

90° 

2000  m. 

583,922 

14.2 

90° 

500  m. 

1,397,848 

21.9 

Azimuth  and  Coordinate  Results 


Azimuth 

90° 

Coordinate 

X 

4 

0.1 

60° 

X 

555 

0.4 

45° 

Z 

10,838 

1.9 

90° 

Z 

39,780 

3.7 

60° 

Y 

71,866 

5.0 

60° 

Z 

120,449 

6.4 

45° 

X 

173,386 

7.7 

45° 

Y 

205,921 

8.4 

90° 

Y 

3,165,470 

33.0 
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TABLE  10  -  (CONT'D) 


Azimuth  and  Resolution  Results 
Azimuth  Resolution 


1 

45° 

10  m. 

974 

0.4 

2 

90° 

5  m. 

5576 

0.9 

3 

60° 

5  m. 

8948 

1.2 

4 

45° 

5  m. 

12159 

1.4 

5 

60° 

10  m. 

246,983 

6.1 

6 

90° 

10  m. 

2,271,297 

18.6 

Azimuth  and  Range  Results 
Azimuth  Range 


1 

90° 

50  km 

412 

0.4 

2 

450 

100  km 

1,153 

0.6 

3 

60° 

100  km 

2,203 

0.9 

4 

60° 

50  km 

6,767 

1.5 

5 

45° 

50  km 

19,367 

2.6 

6 

90° 

100  km 

23,482 

2.8 

7 

45° 

200  km 

98,991 

5.8 

8 

60° 

200  km 

392,999 

11.6 

9 

90° 

200  km 

1,983,070 

26.1 

Range  and  Resolution  Results 
Range  Resolution 


1 

100  km 

5  m. 

93,358 

3.8 

2 

50  km 

10  m. 

132,844 

4.5 

3 

100  km 

10  m. 

222,553 

5.8 

4 

50  km 

5  m. 

319,910 

7.0 

5 

200  km 

5  m. 

1,324,096 

14.2 

6 

200  km 

10  m. 

1,437,712 

14.8 

TABLE  11 

-  AVERAGE  ERRORS 

FOR  INTERACTIONS 

BETWEEN 

THREE  OF  THE 

MAIN  FACTORS 

Order  of 

Results 

Combination  Values 

Ave. 

Error  (-  m. ) 

Range  (km) 

Terrain  Ht  (m) 

Coord. 

1 

50 

2000 

X 

0.4 

2 

50 

500 

X 

1.1 

3 

50 

0 

X 

1.8 

4 

50 

0 

Y 

2.0 

5 

100 

500 

X 

2.0 

6 

200 

500 

Y 

2.1 

7 

50 

0 

Z 

2.1 

8 

100 

2000 

Y 

2.2 

9 

100 

2000 

X 

2.5 

10 

50 

500 

Z 

2.5 

11 

50 

500 

Y 

3.5 

12 

200 

500 

Z 

4.0 

13 

50 

2000 

Y 

4.1 

14 

100 

2000 

Z 

6.1 

15 

100 

0 

Z 

8.1 

16 

100 

0 

X 

8.2 

17 

100 

0 

Y 

9.6 

18 

100 

500 

Z 

9.6 

19 

100 

500 

Y 

10.2 

20 

50 

2000 
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3.5  DISCUSSION  OF  RESULTS 

The  purpose  of  the  study  was  to  determine  the  relative  significance 
of  several  mapping  variables  used  in  locating  terrain  features  with 
radar  and  DTM  data.  Assumed  standard  errors  were  employed  which  may  or 
may  not  be  applicable  for  specific  cases.  In  the  long  run  these  specific 
error  magnitudes  should  not  affect  the  relative  significance  of  the 
factors  involved.  In  addition  the  relative  sizes  of  the  errors  in  Tables 
9,  10  and  11  should  not  be  greatly  changed  provided  the  errors  generated 
are  randomly  distributed. 

Considering  the  five  main  factors  it  was  shown  that  the  coordinate 
factor  was  significant.  The  actual  average  errors  for  the  21 -point 
scenes  provided  the  input  for  the  ANOVA  computation.  Based  on  this 
computation  there  is  a  standard  deviation  of  121  m.  which  can  be  attri¬ 
buted  to  the  variation  between  X,  Y  and  Z  coordinates.  The  standard 
deviation  attributed  to  the  experimental  error  (repeating  the  computation 
three  times  with  different  random  errors)  was  75  m.  The  weakness  in 
consistency  for  coordinate  computation  is  largely  attributed  to  poor 
geometry  and  weak  functional  relationships  during  the  computations.  A 
comparison  of  Tables  3  and  4  shows  an  example  where  an  azimuth  of  45° 

(Table  3)  provides  approximately  equal  standard  deviations  for  X  and  Y 
but,  where  the  azimuth  changes  to  90°  (Table  4),  the  standard  deviation  for 
Y  is  approximately  three  times  that  for  X.  The  values  for  average  error 
are  also  greater  for  the  Y  coordinate  as  seen  in  Table  9. 

Although  range  shows  similar  tendencies  the  range  standard  deviation 
of  108  m.  was  not  significant  at  the  10%  level.  Comparing  Tables  2  and 
3  it  can  be  seen  that  the  standard  deviations  for  X,  Y  and  Z  coordinates 
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are  several  times  larger  at  a  range  of  200  km.  as  opposed  to  50  km. 

In  addition.  Table  9  shows  the  overall  average  error  to  be  much  greater 
for  the  200  km.  range. 

The  standard  deviation  for  resolution  is  80  m.  which  was  not  signifi¬ 
cant  when  compared  to  the  experimental  error  of  75  m.  The  effect  of 
resolution  can  be  seen  by  comparing  Tables  3  and  6.  The  average  errors 
and  standard  deviations  are  generally  larger  for  the  10  m.  resolution. 

Terrain  height  variations  yield  a  standard  deviation  of  75  m.  which 
is  not  significant.  Tables  3  and  5  show  that  a  change  from  0  to  2000  m. 
for  average  terrain  height  will  increase  the  average  errors  and  standard 
deviations  for  all  coordinates. 

The  least  significant  of  the  main  factors  was  azimuth  whose  varia¬ 
tion  accounted  for  a  standard  deviation  of  60  m.  Tables  3  and  4 
show  the  error  pattern  is  more  consistent  and  generally  smaller  for  an 
azimuth  of  45°  as  compared  to  the  90°  azimuth. 

The  five  main  factors  can  be  combined  into  10  different  pairs. 

Table  8  indicates  that  the  pairing  of  range  and  terrain  height  is  signif¬ 
icant  with  a  standard  deviation  of  110  m.  Some  of  the  9  combinations  of 
range  and  terrain  height  covered  in  this  study  are  to  be  avoided.  Table 
10  indicates  that  the  combination  of  200  km.  range  and  2000  m.  terrain 
height  is  the  main  contributer  to  the  significance  of  this  pairing. 

Refer  to  Figure  7  for  an  overview  of  the  average  errors  produced  by  the 
range-terrain  height  combinations.  Each  of  the  9  plotted  points  is  an 
average  error  compiled  from  54  values  representing  1134  rays. 

When  considering  the  combination  of  range  and  coordinate,  the 
standard  deviation  due  to  the  9  possible  combinations  was  93  m.  This 
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combination  was  not  significant.  Table  10  reveals  that  most  of  the  vari¬ 
ation  is  attributed  to  the  Y  coordinate  at  200  km.  Figure  8  displays  the 
coordinate  trends  at  the  three  ranges.  Once  again  54  values  are  repres¬ 
ented  by  each  plotted  point.  The  trends  in  X  and  Y  are  fairly  consistent 
but  Z  displays  some  randomness. 

The  9  pairings  of  terrain  height  and  coordinates  were  again  not 
significant  with  a  standard  deviation  of  91  m.  Table  10  gives  the  results 
of  the  pairings  which  are  plotted  in  Figure  9.  The  most  consistent  results 
for  all  coordinates  were  obtained  at  the  500  m.  terrain  height.  The  X 
coordinate  was  the  most  consistent. 

The  variations  ascribed  to  azimuth  and  terrain  height  combinations 
were  basically  the  same  magnitude  as  those  introduced  by  the  experimental 
error.  Table  10  lists  the  average  errors  which  are  plotted  in  Figure 
10.  Fairly  consistent  trends  are  exhibited  and  low  terrain  heights  and 
45°  to  60°  azimuths  are  favored. 

The  results  for  three  factor  combinations  showed  that  range,  terrain 
height  and  coordinates  were  a  significant  combination  which  produced  a 
standard  deviation  of  104  in.  between  the  27  combinations  possible.  Table 
11  lists  the  average  errors  which  each  represent  18  values.  It  can  be 
seen  that  range  values  behave  fairly  consistently  with  best  results  at 
50  km.  and  poorest  results  at  200  km.  The  terrain  heights  results  are 
fairly  random  throughout  the  list  of  27  and  the  coordinate  listings 
generally  favor  X,  Y  and  Z  in  that  order.  Most  other  three  factor 
combinations  produce  less  variation  than  the  experimental  error  and 
consequently  are  not  at  all  significant. 

One  of  the  four  factor  combinations  was  significant,  namely; 
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Terrain  Height  (m.) 


-  AVERAGE  ERROR  VALUES  FOR  TERRAIN  HEIGHT  AND 
COORDINATE  COMBINATIONS. 


Terrain  Height  (m. ) 


-  AVERAGE  ERROR  VALUES  FOR  TERRAIN  HEIGHT 
AND  AZIMUTH. 


resolution,  terrain  height,  coordinate  and  azimuth  (CDFA).  All  other 
four  factor  combinations  did  not  approach  significance.  Fifty-four 
values  are  represented  in  the  CDEA  combination.  Detailed  results  show 
the  greatest  variation  is  caused  by  10  m.  resolution,  500  m.  terrain 
height,  Z  coordinate  and  90°  azimuth;  however,  only  nine  values  or  189 
individual  rays  are  represented  by  this  result. 

The  particular  significance  of  the  results  will  depend  on  the  appl 
cation,  however,  specific  data  collection  designs  are  necessary  to  fine 
tune  the  mapping  system.  For  this  particular  set  of  data  the  results 
indicate  that  one  should  collect  radar  data  at  an  azimuth  of  45°,  a 
range  of  100-50  km,  a  resolution  of  5  m.  and  over  low  level  terrain  for 
the  best  results. 

4-0  CONCLUSIONS  AND  RECOMMENDATION 
4.1  CONCLUSIONS 

The  following  conclusions  are  based  on  the  fixed  factors  of  this 
experiment. 

1.  Avoid  collecting  radar  data  when  the  axes  of  the  radar  coordi¬ 
nate  system  are  parallel  with  axses  of  the  DTM  coordinate  system 
DTM  coordinates  along  track  (Y)  contained  large  errors  because 
of  this.  If  unavoidable,  an  intermediate  axes  rotation  may 
prove  helpful. 

2.  Radar  operating  at  maximum  range  generally  yields  less  reliable 
measurement  results. 

3.  A  small  increase  in  resolution  greatly  improves  the  geometric 
integrity  of  the  mapping  effort. 
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4.  Increased  terrain  height  leads  to  decreased  accuracy  and  pre¬ 
cision  in  position  determination. 

5.  Variable  geometry  results  in  considerable  variation  in  the  relia¬ 
bility  of  X,  Y,  and  Z  coordinate  determination.  The  conditions 
in  this  experiment  were  most  unfavorable  for  the  Y  coordinate. 

6.  The  proper  design  of  radar  flights  and  data  reduction  is  essen¬ 
tial  if  the  data  is  to  be  used  to  favorably  compliment  mapping 
efforts. 

4.2  RECOMMENDATION 

The  experiment  identified  important  parameter  variations  and 
dependencies.  More  extensive  analysis  of  data  collection  parameters 
would  prove  productive;  particularly  when  sensors  are  combined  for  a 
common  cause.  Computer  simulations  are  relatively  inexpensive  yet  often 
pin  point  the  not  so  obvious.  It  is  recognized  that  simulations  are  far 
removed  from  real  world  data  collection,  analysis  and  mapping;  however, 
experimental  data  is  available  for  establishing  variance  on  measurement 
systems  and  this  data  can  be  used  to  fortify  simulation  results. 

The  combination  of  simulation  and  factorial  analysis  has  been  ap¬ 
plied  successfully,  and  it  is  recommended  that  this  data  reduction  com¬ 
bination  be  used  more  frequently  in  the  preliminary  design  of  remote 
sensing  systems. 
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APPENDIX  A 


BASIC  PROJECTION  EQUATIONS  FOR  REAL-APERTURE  RADAR 


The  purpose  of  this  appendix  is  to  explicitly  set  down  the  basic 
projection  equations  for  real  aperture  radar  as  used  in  this  error 


propagation  study.  Although  these  equations  have  been  published  in 
several  forms  (Leberl,  1972  and  Leber! ,  1976  particularly)  their  avail¬ 
ability  is  still  limited  and  it  is  felt  that  a  development  of  the 
versions  used  in  this  report  would  be  helpful. 


Real -aperture  radar  imagery  can  be  expressed  in  terms  of  a  gen¬ 


eral  projection  equation  as  given  in  A-l,  A-2,  and  A-3  below. 


(See  Figure  A-l ) . 


Xp  =  XR  +  AT  BT  •  xp  +  AX  (A-l) 
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FIGURE  A-l  COORDINATE  SYSTEM  FOR  RADAR  PROJECTION  EQUATIONS 


A  simplification  of  equation  (A-3)  will  now  be  made  by  assuming 
that  the  sensor  is  parallel  and  coincident  with  the  mount.  In  this 


case  B  becomes  a  unit  matrix  and  aK  is  zero. 

xp 

yp 

ZP 

where:  xp  =  (sin  t)  r 

yp  =  (sin^  n-sin^-r)'5  r 
zp  =  (-  cos  ft)  r 

r  =  Range  from  radar  to  DtM  point  P. 
t  =  Squint  angle 
ft  =  Elevation  angle 

The  A  matrix  utilized  in  this  Study  in  (A-4)  is  defined  below: 

cos  <p  cos  k  ;  cos  w  sin  k  +  sin  w  sin  <j>  cos  k  ; 

cos  <p  sin  k  ,  cos  w  cos  k  -  sin  w  sin  <j>  sin  k  ' 

sin  <p  ;  -  sin  w  cos  t  1 

sin  w  sin  k  -  cos  w  sin  <t>  cos  k 
sin  w  cos  k  +  cos  w  sin  <)>  sin  k 
cos  <p  cos  w  _  (A-5) 

+w  =  clockwise  rotation  about  x  axis  -  primary  rotation 
+<)>  =  "  "  "  y  axis  -  secondary  rotation 

+k  =  "  "  "  z  axis  -  tertiary  rotation 

(Figure 


where: 


In  order  to  provide  a  complete  sample  computation  such  as  that 
carried  out  in  the  study  it  will  be  necessary  to  refer  to  a  sketch  of 
the  ficticious  locations  of  radar  and  ground  position.  (See  Figure  A-2). 

R  represents  the  radar  location  in  an  aircraft.  The  coordinates 
of  the  radar  in  the  DTM  system  are: 

XR  =  0  m.;  Yr  =  +  99,254  m. ;  ZR  =  +  12,192  m. 

P  represents  the  target  or  point  of  interest  in  the  DTM  system. 

The  coordinates  of  P  are: 

Xp  =  0  m. ;  Yp  =  0  m. ;  Zp  =  0  m. 


Utilizing  these  coordinates  range  is  computed  to  be  100,000  m. 
The  exterior  orientation  of  the  radar  is: 


K  =  0  degrees;  cj>  =  0  degrees;  w  =  0  degrees 

(A-4)  can  be  rearranged  as  (A-6)  and  utilized  to  compute  t  (squint 
angle)  and  n  (the  elevation  angle). 
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Substituting  the  values  for  our  example  yields: 


(sin  t)  100,000 

1  0  0 

0  -  0 

"o  i 

(sin2n  -  sin2  x)  \  100,000 

= 

0  1  0 

0  -  99,254 

= 

-99,254 

(-  cos  a)  100,000  _ 

0  0  1 _ 

0  -  12,192 

-12,192 

yP 


Solving  for  t  and  fi: 

0 


t  =  sin-1 


=  cos 


-1 


100,000 

-12,192 


=  0  degrees 


=  83  degrees 


L  -100,000 

This  completes  the  computation  of  a  single  ficticious  ray  from 
radar  to  ground  point.  Each  ray  was  then  perturbed  with  random  errors. 
Errors  were  introduced  for  radar  position  (XR,  YR,  ZR),  radar  orientation 
(k,  4» ,  w),  and  ray  length  and  direction  (r,  t,  n). 

For  example  assume  the  following  random  errors  occurred  in  the 
above  9  parameters: 


R 


=  +  9  m. 

=  -3  m. 


ZR  =  +20  m. 


k  =  +0.1  deg. 
<P  =  -0.1  deg. 

w  =  -0.1  deg. 


r  =  +  5  n.. 
t  =  +  0.003  deg. 
a  =  -  0.003  deg. 


These  error  magnitudes  were  considered  reasonable  values  from 
previous  studies.  The  adjusted  parameter  values  are  now: 
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9  m. 

k  =  +0.1  deg. 

r 

=  100,005  m. 

99,251  m. 

4>  =  -  0.1  deg. 

T 

=  +0.003  deg. 

12,212  m. 

w  =  -  0.1  deg. 

a 

=  82.997  deg. 

Utilizing  equation  (A-4)  the  coordinates  of  P  now  are: 

Xp  =  -  137.72  m.i  Yp  =  -28.91  m. ;  Zp  =  192.85  m. 

Recall  the  correct  value  for  each  coordinate  was  zero. 

The  total  experiment  entailed  the  computation  and  perturbation 
of  ficticious  rays  or  162  scenes  each  having  21  rays  or  points  of  inter¬ 
est. 

These  coordinate  errors  were  the  bases  of  the  input  for  the  analysis 
of  the  variance  computation. 
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APPENDIX  B 


ANALYSIS  OF  VARIANCE  MODEL 


The  factors  studied  in  this  experiment  were  azimuth,  range, 
resolution,  terrain  height,  and  coordinates.  The  purpose  of  the 
analysis  of  variance  procedure  was  to  aid  in  selecting  the  best 
combinations  of  the  above  variables  when  using  radar  to  locate  points 
of  interest  on  the  terrain.  Assuming  that  all  factors  in  the  error 
propagations  are  constant  except  those  above,  there  is  a  need  to  analyze 
the  propagated  errors  in  such  a  way  that  one  can  determine  not  only 
how  each  individual  factor  affects  the  results  but  in  addition,  how 
the  interactions  of  all  factors  affect  the  results.  A  factorial  exper¬ 
iment  was  used  to  accomplish  this.  The  levels  of  each  factor  were 
selected  to  introduce  ranges  in  the  variables  which  were  of  particular 
interest  in  this  study.  This  produced  a  fixed-effects  model  which  is 
shown  schematically  in  Figure  B-l. 

The  mathematical  model  which  represents  the  factorial  experiment 


is: 
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ijklmno  =  propagated  error 

A  =  azimuth-fixed  D  =  terrain  height-fixed 

B  =  range-fixed  E  =  coordinate-fixed 

C  =  resolution-fixed  6  =  experimental  error 

a . ^  is  the  interaction  effect  between  the  ith  azimuth  and  the  jth 

range,  etc.  The  number  of  levels  for  each  of  the  fixed  factors  are: 

I  =  3,  J  =  3,  K  =  2,  L  =  3,  and  M  =  3. 

It  is  assumed  that  the  linear  additive  model  as  written  above  is 

A 

applicable  for  this  study.  The  a‘s  are  fixed  effects  and  Ea.  = 

£«j  =  "  0  is  assumed  true'  The  61jklmnos  are  the 

experimental  unit  and  they  are  assumed  to  be  normally  and  independently 

2 

distributed  with  zero  mean  and  variance  . 

For  each  item  in  the  model  equation,  the  following  type  of  hypoth¬ 
esis  is  tested. 
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REPETITION 


FIGURE  B-l  ANALYSIS  OF  VARIANCE  TABLE  FOR  MAIN  EXPERIMENT 
EACH  CELL  CONTAINS  AVERAGE  ERROR  FOR  21  TERRAIN  POINTS. 
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